Periodic Landau gauge and Quantum Hall effect in twisted bilayer graphene 



m 

o 

(N 

Oh! 

< 

(N 



Yasumasa Hasegawa^ and Mahito Kohmoto^ 

^Department of Material Science, Graduate School of Material Science, University of Hyogo, 

3-2-1 Kouto, Kamigori, Hyogo, 678-1297, Japan 

^Institute for Solid State Physics, University of Tokyo, 

5-1-5 Kashiwanoha, Kashiwa, Chiha 277-8581, Japan 

(Dated: April 29, 2013) 

Energy versus magnetic field (Hofstadter butterfly diagram) in twisted bilayer graphene is studied 
theoretically. If we take the usual Landau gauge, we cannot take a finite periodicity even when the 
magnetic flux through a supercell is a rational number. We show that the periodic Landau gauge, 
which has the periodicity in one direction, makes it possible to obtain the Hofstadter butterfly 
diagram. Since a supercell can be large, magnetic flux through a supercell normalized by the 
flux quantum can be a fractional number with a small denominator, even when a magnetic fleld 
is not extremely strong. As a result, quantized Hall conductance can be a solution of nontrivial 
Diophantine equation. 

PACS numbers: 73.22.Pr, 73.20.-r, 73.40.-c, 81.05.ue 
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I. INTRODUCTION 

Two dimensional electron systems are realized in 
grapheneii Bilayer graphene^ and twisted bilayer 
graphene^i have been shown to have many interesting 
properties and been studied extensively. The quan- 
tum Hall effects in single layer graphene^ and bilayer 
graphene are well understood by the energy versus mag- 
netic field, which is known as the Hofstadter butter- 
fly diagram, in the tight binding models of honeycomb 
latticefir^, and on Bernal-stacked bilayer graphene^.. The 
quantized value of the Hall conductance is obtained by 
the solution of the Diophantine equatiouiiSrJ^ Near half 
filled case, the quantum Hall conductance is given by a 
solution of the Diophantine equation Eq. ([24| with s,. = 1 



rrij = 2, m^ = 1, L^ = m^a^ + m^ a^ 






f^Ti / 3 - a 




FIG. 1. (color online). Twisted bilayer graphene with 
(mi, 7712) — (2, 1). Large (small) open and closed circles are A 
and B sublattice in the first (second) layer, respectively. The 
second layer is rotated by 7r/3 — a. The orange area is the 
supercell. 



{sr = 2) in single layer (bilayer) graphene. 

Twisted bilayer graphene attracts much attention re- 
cently. When two layers are twisted in a commensurate 
way, a supercell becomes large with a moire pattern^ii^"— 
(see Fig. [ij and the velocity at the Dirac points is shown 
to become smaller when the rotation angle a is smal h^^i^^ 

Quantum Hall effect and the Hofstadter butterfly di- 
agram in moire superlattices have been observed exper- 
imentally in single layer graphene on hexagonal boron 
nitride (hBN)^^ and Bernal-stacked bilayer graphene on 
hBN.2" 

Electronic structure in twisted bilayer graphene in 
a uniform magnetic field has been studied only by 
taking the linearized energy dispersion near the Dirac 
points^"— or by the Lanczos algorithm applied to large 
systems in real space^i. The whole lattice structure has 
not been taken into account when we take the linearized 
energy dispersion. Exact band structure is difficult to 
obtain by the Lanczos algorithm due to the finite-size 
effects and numerical errors. 

When we use the usual Landau gauge in twisted bilayer 
graphene with long range hoppings, there are no period- 
icity in the phase factor of hoppings. In that case we 
cannot obtain the Hofstadter butterfly diagram. In this 
paper we show that we can recover the periodicity when 
a magnetic flux through a supercell is a rational number, 
if we use the periodic Landau gauge. We calculate the en- 
ergy spectrum in bilayer graphene in a magnetic fleld in a 
similar way as in single layer graphene or Bernal stacked 
bilayer graphene. We obtain very rich Hofstadter dia- 
grams. There appear many energy gaps, which can be 
labeled by integers of solutions of the Diophantine equa- 
tion. The quantized value of the Hall conductance is 
obtained by these integers. 



II. TWISTED BILAYER GRAPHENE 

In a unit cell in each layer there are two sites, A and 
B, which form triangular lattices respectively. We define 
unit vectors as 



SLi = a 



v2 
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_1 
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and 



a2 — Rtt/3^1 



%/3 



(1) 



(2) 



where a is the lattice constant and Rtt/^sli is the 7r/3 ro- 
tated vector of ai . Hereafter we take a — 1 for simplicity. 
The reciprocal lattice vectors are given by 
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In the first layer, sites in the A sublattice are given by 
sets of two integers (ji,J2) as 



Jiai + J2a2. 



(5) 



Three vectors connecting nearest neighbor sites in the 
first layer are 



Sa = 
Sb = 
Sr = 



V3 
3 




V3 
,6 



V3 



(6) 

(7) 
(8) 



Sites in the B sublattice in the first layer are given by 






Sa. 



(9) 



The AB (Bernal) stacking of bilayer graphene is obtained 
by rotating the second layer around one of the A site 
in the first layer by the angle (2n + l)7r/3, where n is 
an integer. In this case the A sublattice in the second 
layer is just above the A sublattice in the first layer, 
but the B sublattice in the second layer is on the cen- 
ter of the hexagon in the first layer. The same stacking 
is obtained by translating the second layer by —Sa^ —Sb 
or —Sc- When the rotation angle is 2mr/3, we obtain 
the AA stacking, i.e., all sites in the second layer are 
on the sites in the first layer. We obtain twisted bilayer 
graphene, when the rotation angle is neither {2n+ l)7r/3 
nor 2n7r/3. 

When twisted bilayer graphene has supercell with finite 
number of sites, it is called commensurate twisted bilayer 
graphene. We construct commensurate twisted bilayer 



graphene as follows; Since there is six-fold symmetry in 
twisted bilayer graphene, we can take a supercell as a 
diamond with the angle 7r/3 as shown in Fig. [TJ We 
define unit vectors of superlattice with two integers toi 
and 77i2 {rni ^ 0, TO2 ^ 0, and \mi\ ^ |'Ti2|): 



T 1 = r-^ 

1 7711,7712 

L2 = Rtt/3^1 



rriiai 



™2a2, 

■7712, (7711+7712)' 



(10) 

(11) 



Twisted bilayer graphene with (mi , 7712) 
in Fig. [TJ Since 



ai • 3.2 



(2,1) is shown 



(12) 



we obtain 

|LiP = |L2| 
The area of a supercell is given by 



nil + "''2 + T?,iTO2 = Uq. 



S 



ILillL, 



V3 



-no- 



(13) 



(14) 



There is another site in the supercell that has the same 
distance from the origin as |Li|, which we define L[ as 
shown in Fig. [1] 



L' 



TO2ai -|-mia2. 



(15) 



We define a by the angle between the vectors h'^ and Li. 
Since 



L'j • Li = -{ml +ml + 2miTO2), 



we obtain 



cos a = 



li'i ■ Li ml + 4miTO2 + m^ 



im|Li| 2{ml+mim2 



(16) 



(17) 



Then we obtain twisted bilayer graphene by rotating the 
second layer with the angle 7r/3 — a to move the vector 
h'l into L2. We obtain another type of twisted bilayer 
graphene when we rotate the second layer by the angle 
—a. In this paper we take the rotation angle 7r/3 — a to 
obtain the Bernal stacking when a — >■ 0. 

The Bravais lattice of twisted bilayer graphene is the 
— a/2-tilted two-dimensional triangular lattice with the 
primitive vectors Li and L2. A supercell has uq sites in 
A and B sublattice in each layer, hence Auq sites. 



III. TIGHT BINDING MODEL AND PERIODIC 
LANDAU GAUGE 



We consider tight binding models in a uniform mag- 
netic field. Spin is not taken into account. When a mag- 
netic field is applied, the hopping t^ between sites r^ and 
Tj [vi and Ti are on the same layer or on different layers) 
has the factor exp(z0y) with a phase 9ij given by 



2n 



A-dL 



(18) 



where A is a vector potential and 

00 = — , 
e 



single layer graphene (nearest-neighbor hoppings) 



(19) 



is the flux quanta with charge e, the speed of light c and 
the Planck constant h. The Hamiltonian is 



n 






^ijC^j^i 



h.c. 



(20) 



J 



where c] and Ci are the creation 
tion operators at site i, respectively, 
approximatio n-'^^'^'^i^^ , 



and annihila- 
We take the 



t exp 



Hj — < 



ti2 exp - 



d-\S^ 



if sites i and j are on the same layer 
d — cIq^ 



312 



where d 



if sites i and j are on different layers 

(21) 
Tjl is the distance between sites i and 



j and do is the distance between layers. When we take 
(5 — >■ and ii2 — 0, we obtain two independent layers of 
honeycomb lattice with only nearest-neighbor hoppings. 
When ti2, S and (5i2 are finite, we obtain twisted bilayer 
graphene with finite range hoppings. Interlayer hoppings 
are not restricted to the perpendicular direction. 

Note that if we change the signs of cj and Ci in the 
second layer and the sign of ii2 simultaneously, we obtain 
the same Hamiltonian. This means that the energy is 
independent of the sign of ti2. 

Even if the flux per supercell is an integer times the 
flux quantum 4>q, the phase factor 0ij is not periodic 
with modulus 27r, if we use the usual Landau gauge 
(A = Hxy). For single layer graphene with only nearest- 
neighbor hoppings, we could take a special gang' 



6.7.26 



which the phase factor appears only in the links for one of 
the three directions, Sa, Sb or Sc- However, such choice of 
gauge is not possible for twisted bilayer graphene. In or- 
der to obtain the energy spectrum of tight-binding mod- 
els on twisted bilayer graphene in a magnetic field, we 
have to choose a suitable gauge.^^^"— ° As far as we know, 
the special choice of the gauge was first used in the system 
of the 4x4 square lattice with periodic boundary condi- 
tions in the presence of the uniform magnetic flux p/16 
with p = 1,2,3,---.— We explain the periodic Landau 
gauge for the square lattice in Appendix [Xj The general- 
ization to the non-square lattice is given in Appendix IBJ 
When flux through a supercell is 



9 



(22) 



where p and q are mutually prime integers in twisted 
bilayer graphene with commensurate twisted angle 
fEg. dTTl) '). energy spectrum is obtained by the eigenval- 
ues of (4nog) X (Anoq) matrix, as in the case of the single 
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FIG. 2. (color online). Energy spectrum in single layer 
graphene with only nearest-neighbor hoppings. Numbers in 
the figures are {sr,tr). Quantized value of Hall conductance 
is given by U. See Eqs. (pi)) and (f25l) . 



layer graphene where it is obtained by the eigenvalues of 
(2g) X {2q) matrix.^ 

In Figs. [21 we plot the energy versus magnetic flux 
through a unit cell in single layer graphene with only 
nearest-neighbor hoppings. In Figs. [3l HI and [5l we 
take parameters for the bilayer graphene^^ t = 2.7 eV, 
ti2 = -0.48 eV (tiz/t = -0.18), l^^l = 0.142 nm, do = 
0.335 nm, and S ^ S12 ^ 0.0453 nm {6/\Sa\ = 0.184), and 
we plot the energy versus magnetic flux through a unit 
cell in each layer {(j) = $/no) in twisted bilayer graphene. 



IV. DIOPHANTINE EQUATION 

Consider the case 



$ p 

00 g' 



(23) 



m =3 m=2 n = 19 a =13. 17 q=17 



m^ = 8 m^ = 7 iig = 169 a =4.41 q=5 
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FIG. 3. (color online). Energy spectrum in twisted bilayer 
graphene with (mi, 7712) = (3,2). Numbers in tlie figures are 
(sr,tr)- Quantized value of Hall conductance is given by tr- 



FIG. 4. (color online). Energy spectrum in twisted bilayer 
graphene with (7711,7712) = (8, 7). Numbers in the figures are 
(sr,ir)- Quantized value of Hall conductance is given by tr- 



where $ is the flux through a supercell, p and q are in- 
tegers. In square lattice and honeycomb lattice $ = 0, 
and in twisted bilayer graphene $ = no^, where (/) is the 
flux through a unit cell in each layer. When chemical 
potential is in the rth gap from the bottom, we have the 
Diophantine equationiSr— i^ 



r = qsr 



ptr 



which gives quantized Hall conductance by 






(24) 



(25) 



If we take account of the spin and neglect the Zeeman 
energy, the Hall conductance is multiplied by 2. 

For the tight binding models with only nearest- 
neighbor hoppings in square lattice or honeycomb lattice, 
the flux quantum c/jq through a unit cell is equivalent to 
zero magnetic flux. As a result, the energy spectrum 
is periodic with respect to (f) with a period (f>o- Even 



when we consider the models with long range hoppings, 
the energy spectra are periodic function of </> with a pe- 
riod 200 or 600 in the square lattice or the honeycomb 
lattice, respectively. This is because the smallest areas 
enclosed by hoppings are 1/2 and 1/6 of the areas of a 
unit cell in the square lattice and the honeycomb lattice, 
respectively. See Fig. [5] for the honeycomb lattice. The 
energy spectrum is also periodic with respect to (j) with a 
period 60o for Bernal stacked bilayer graphene. The sit- 
uation is drastically changed in twisted bilayer graphene. 
When there are hoppings between layers in twisted bi- 
layer graphene, projected areas enclosed by hoppings 
have irrational values as shown in the red triangles in 
Fig. [6l As a result the energies are not periodic in 0. 

In single layer graphene, there are 2q band when flux 
through a supercell is {p/q)(j)o- In Fig. [51 we show {sr,tr) 
for several gaps for single layer graphene^. Large gaps 
have indices tr = 0, 1 , and 2. The gaps, which are 
focused at the bottom of the band at — >■ 0, have Sr = 
0, and tr = 1,2,3, •••. They correspond to the usual 
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FIG. 5. (color online). Energy spectrum in twisted bilayer 
graphene with (mi,m,2) = (12,11). Numbers in the figures 
are (sr,ir)- Quantized value of Hall conductance is given by 



Landau levels. The gaps, which are focused at the top of 
the band at (^ — )■ 0, have Sr = 2, and tr = —1, —2, —3, • • • . 
The gaps near half filling [E sa 0) and <C 0o have Sr = 
1 and tr = ±1,±3,±5,---, which have been observed 
in graphene^i^. For the finite (f) the band near E k, Q 
becomes broadened gradually and many gaps can be seen 
in Fig. [21 Note that = 0o in a unit cell corresponds 
to 40, 000 T, which is not attainable in a present day 
laboratory. 

In Figs. [21 [H and [SI we plot the Hofstadter butterfly 
diagrams for twisted bilayer graphene with (7711,7712) = 
(3,2), (8,7) and (12,11), respectively. Although energy 
gaps near the bottom have Sr — and i^ = 1, 2, 3, • • • , for 
all three cases as in single layer graphene, there are cross- 
ings of the bands near the bottom of the energy. For ex- 
ample, the gap indexed by (0, 2) vanishes at (/)/0o ~ 0.16 
and E/t « —2.8, at which band-crossing occurs. These 
crossings of bands can be understood by the indepen- 
dent Landau levels for the two local minimums of the 
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FIG. 6. (color online). Twisted bilayer graphene. When there 
are only nearest-neighbor hoppings in each layer, an area en- 
closed by hoppings is S (blue hexagon). When there are 
second-nearest-neighbor or third-nearest-neighbor hoppings 
in each layer, areas enclosed by hoppings are S/6 and its mul- 
tiples (blue triangle). The hoppings between layers make the 
area enclosed by hoppings to be irrational number times S 
(red triangles). 





FIG. 7. (color online). Schematic band structure near the 
bottom of the energy. There are two local minimums and 
nearly independent Landau Levels (dashed blue lines and solid 
red lines). 



energy in the absence of a magnetic field (see Fig. [7]). 
Near the top of the energy the large gaps are indexed as 
{sr,tr) = (4no,-2),(4no,-4),(4no,-6),-- •, which can 
be understood by the fact that there are 4no bands and 
nearly degenerate two local maxima of the energy in the 
absence of a magnetic field. 

A very interesting feature is seen near half filling. Be- 



sides the large gaps of (s,-,ir) = (2no,±2), many new 
gaps become visible as a becomes small. For example, 
(s,,i,) = (2no±2,0),(2no±2,±2),(2no±4,±2),(2no± 
6, ±2), • • • are seen in the lower figures in Figs. S] and [S] 
These new gaps are caused by a large supercell, which 
has 4no sites. Since Hall conductance is given by i^, 
the band between the gaps with same t^ ((2no,2) and 
(2no — 2, 2), for example) does not contribute to the Hall 
conductance. Mathematically, that band is a Cantor set 
and consists of narrower bands and much smaller gaps. 
Each narrow band gives a finite Hall conductance and 
the total contribution vanishes. 



V. SUMMARY 

We obtain the Hofstadter butterfly diagram for twisted 
bilayer graphene. The use of the periodic Landau gauge 
is crucial. Due to large number of sites (4no) in a su- 
percell, a rich structure of the Hofstadter butterfly di- 
agram appears, especially near half- filling and for small 
rotation angle a. The gaps are indexed by two integers 
Sr and tr (Eq. ([M)) ). The Hall conductance is given 
by Oxy — {e^ /h)tr. While gaps with s^. = 0, 1 and 
2 are large in single layer graphene, many gaps with 
Sr = 2no, 2no ± 2, 2no ± 4, • • • become large as a be- 
comes small. Since a supercell becomes large as a — >■ 0, 
flux per supercell versus the flux quantum ($/(?!)o — p/q) 
can be a rational number with a small denominator q in 
not an extremely strong magnetic field. For mi =^ 12 
and 7712 = 11, we obtain 77o = 392. In that case 100 
T corresponds to ^/4>q ~ 1, which may be attained in 
experiment. 

Near half-filling, there are many narrow bands, which 
do not contribute to the Hall conductance when they 
are completely filled. The narrow bands, in fact, consist 
of even narrower bands, since an energy spectrum is a 
Cantor set. Therefore one can expect non-trivial Sr and 
tr in the narrower band. It may be possible to observe 
these phenomena experimentally. 



Appendix A: periodic Landau gauge in square lattice 

We explain the periodic Landau gauge in the square 
lattice. A vector potential A gives a magnetic field 



H = V X A. 



(Al) 



For a uniform magnetic field H = Hi, one can take the 
Landau gauge 



A(^) = Hxy, 



(A2) 



where y is the unit vector along the y direction and there 
is no dependence on y. In this case, however. Ay is not 
periodic in the x direction as shown by the dashed line 
in Fig. [51 We can obtain other vector potential by gauge 
transformation, i.e., adding Vx(r) to A. It is crucial to 
have periodicity in a gauge of twisted bilayer graphene. 
We take 



X(r) = -H\x\y, 



(A3) 



and so, 

A(P^)=i/(xy-V(La;j7j)) 

= H {{x - [x\)y - y ^ (5(a:: - 77 + e)x j , 

\ n— — oo / 

(A4) 

where \_x\ is the floor function (largest integer not greater 
than a;), i.e. x — \_x\ is the fractional part of x. In this 
gauge, which we call the periodic Landau gauge. Ay ' 
is periodic with respect to x with period 1, as shown 
in Fig. [51 Note that Ax depends on y and it is not 
periodic in the y direction. However, the dependence 
of y in A*^^^) appears always with the delta function, so 
exp(70jj) is periodic in the y direction, as we show below. 




FIG. 8. (color online). The y component of vector potentials 
for Landau gauge {Ay , dashed line) and periodic Landau 
gauge {A)f ' , thick blue line) for the square lattice. 



The difference of the usual Landau gauge and the pe- 
riodic Landau gauge is seen in Fig. [S]in which a uniform 
magnetic field with {p/q)4>o through a unit cell is applied 
to square lattice. If we take the usual Landau gauge, 
the phase factor is zero except vertical links as shown in 
Fig- El (a)- The phase factor for the vertical links aX x — n 
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is Bij = l-Knpjq. The periodicity in the x direction is q 
times larger than the periodicity in the absence of a mag- 
netic field. However, if there are other sites in a unit cell, 
the periodicity of the system is changed (this is the case 
in twisted bilayer graphene, where there are 4no sites in 
the supercell (see Figs. [T] and E])). In order to demon- 
strate it in the square lattice, we add sites at (n-l-71, ttt.), 
(n, m -I- 72), and at (n -I- 71, m -)- 72), where n and m are 
integers, < 71 < 1 and < 72 < 1, as shown by the 
filled green circles in Fig.jQjb). The phase factors for the 
links connecting the nearest neighbor sites are shown in 
Fig. Elb) . The phase factor for the vertical link between 
Vi = (n,m) and r^ = {n,m + 72) is 9ij = 2TTnj2p/q- If 
72 is an irrational number, exp(i9ij) cannot be periodic 
with respect to x. 

The periodicity is recovered by taking the periodic 
Landau gauge (Eq. (|A4|) '). The periodicity is 1 in the 
X direction, since A'^^' is periodic in the x direction. 
A uniform magnetic field through each small rectangles 
is ensured by the nonzero phases for the horizontal links 
which come from the delta functions in Eq. (IA4p as shown 
by thick red horizontal arrows in Fig,[SlJc) . The phase fac- 
tor for the horizontal link between r^ = (n-l-71 , to) and 
Tj = {n+1, m) is 9ij = —2Trmp/q, which does not depend 
on 71 . The periodicity in the y direction is q times larger 
than that without magnetic field in the periodic Landau 
gauge. 



Appendix B: periodic Landau gauge in non-square 
lattice 



The periodic Landau gauge discussed in Appendix |X] 
is generalized to non-square two-dimensional lattices, 
which has primitive vectors of the supercell Li and L2, 
which are not orthogonal. The reciprocal lattice vectors 
are 



Fi = 27r 



L2 X z 



(Li X L2) -z' 



and 



2n- 



z X Li 



(Li X L2) -z 
We define oblique coordinate system (^1,^2) by 



(Bl) 



(B2) 



FIG. 9. (color online), (a) The numbers beside the arrows are 
the phases in the usual Landau gauge for the square lattice 
with flux {jp/q)(j>o in a unit cell. The system is periodic in the 
X direction with period q. (b) Sites are added at (n-|-7i,m), 
(71,771-1-72), and (71 + 71,771-1-72), where 71 and m axe integers, 

< 71 < 1 and < 72 < 1 (filled green circles). Blue letters 
are flux in rectangles. The system is periodic only when 72 
is a rational number, (c) The phases (beside arrows) and the 
flux through rectangles (blue letters) in the periodic Landau 
gauge. The system is periodic in the x direction with period 

1 and in the y direction with period q. 



r = 



— ^iLi + C2L2 



(B3) 



where x and y are the coordinates in a orthogonal system. 
For twisted bilayer graphene we have 



Ci = — ^r("^i + 2TO2)a; - niiy 
na \ i 



(B4) 



6 = — ^rl"^! ~ "^2)2; + {mi + TO2)y . (B5) 
no \ 3 / 



The reciprocal vectors are 



and 



_ mi + TO2 , "^2„ 
-T 1 — IJl H *J-2, 

no no 



F2 — Gi H G2. 

no no 



(B6) 



(B7) 



The Landau gauge for the oblique coordinate system is 

^{L,as,) ^ ^^^F2, (B8) 

ZTT 



which is a generalization of Eq. (jA2p to the non-square 
lattice. To make the vector potential periodic with re- 
spect to ^1, we take the periodic Landau gauge as 



A(''^'"^«)=^((a-Laj)F2 

Ztt V 

00 

-6 E '5(6-^ + e)Fi 



(B9) 



In Fig. [To] we show an example of 9ij between the nearest 
neighbor sites in the first layer for mi = 2, m2 ~ 1 
(no = 7) and the flux through a supercell is 0o {p/q — 
1/1). Note that the phase factor 9ij crossing the line 
^1 = n (red numbers in Fig. 1101) is not periodic in the L2 
direction but exp(i6'y) is periodic. 

The advantage of taking the periodic Landau gauge is 
that the vector potential A is periodic in ^1 with period 1. 
This ensures the periodicity of the system in Li direction. 
The system has a periodicity in the L2 direction if the 
flux per supercell is an rational number times ^q. 

If the flux per a supercell is a rational number p/q with 
integers p and q, i.e. 



$ = no0 = ^i/--0o, 

q 



(BIO) 



the phase factor 1 
gauge (Eq. (|B9l) 



kj in Eq. (|18p with the periodic Landau 
has the same value when r^ and r^ are 
translated by Li. It has also the same value when r^ and 
Tj are translated by L2 if the link connecting r.^ and Vj 
does not cross the line ^1 = n, where n is an integer. If 
the link connecting r^ and Vj crosses the line ^1 = n, 9ij 
increases by ±2Trp/q when r^ and Tj are translated by L2 
(± depends on the sign of Fi • d£). Therefore periodicity 
must be q times larger in the L2 direction. 



m^ = 2, m^ = 1 
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FIG. 10. (color online). The first layer of twisted bilayer 
graphene with (mi,m2) = (2, 1) (no — 7). The numbers on 
the links between the nearest neighbor sites show the phase dij 
in the first layer in units of 27r/(18no) = 27r/882 when the fiux 
per supercell is <j>o, i.e. the flux per unit hexagon is 4>o/no (We 
obtain (4- (-35) +424-34- (-155) - 148) x 27r/882 = 27r/7, 
for example). We take the periodic Landau gauge (Eq. (|B9[) V 
The contribution from the delta function in Eq. (|B9|) is finite 
for the red numbers, while it is zero for the blue numbers. 
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